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1. Introduction 

Unbound decaying states are energy eigenstates of a time reversal invariant Hamiltonian 
describing non-dissipative physics in a situation in which there are no particles incident 
PI |2j. Even when the formal Hamiltonian, considered as an operator in the Hilbert 
space of square integrable functions, is Hermitian (self-adjoint), this boundary condition 
makes the eigenvalue problem non- self adjoint and the corresponding energy eigenvalues 
complex, £n = En — zl/2r„ with _E„ > r„ > [Sj. 

Commonly, unbound energy eigenstates are regarded as arising from a perturbation 
with the physics essentially unchanged from the bound state case, except for an 
exponential decay. But, unbound state physics differs radically from bound state physics 
in the presence of degeneracies that is, coalescence of eigenvalues 

In the case of a Hermitian Hamiltonian depending on parameters, the bound state 
energy eigenvalues are real and, when a single parameter is varied, the mixing of two 
levels leads to the well known phenomenon of energy level repulsion [H] and avoided 
level crossing jB]. In the case of unbound energy eigenstates (resonances) of the same 
Hamiltonian, the energy eigenvalues are complex and, when a single parameter is varied, 
this fact opens a rich variety of possibilities, namely, crossings and anticrossings of 
energies and widths [7j, and the so called "change of identity" of the poles of the 
S"— matrix jH]. These novel effects have attracted considerable theoretical [HI HHl HI] 
and recently also experimental interest ^2)- True degeneracies of resonance energy 
eigenvalues result from a joint crossing of energies and widths in a physical system 
depending on only two real parameters and give rise to the occurrence of a double pole 
of the scattering matrix in the complex energy plane PU]- Associated with a double 
pole of the S"— matrix, the quantum system has a chain of Jordan-Gamow generalized 
eigenfunctions ^1 . Examples of double poles in the scattering matrix of simple 
quantum systems have been recently described [13 CEl El HHj • Degeneracy of complex 
energy eigenvalues of non- Hermitian PT symmetric Hamiltonian has been discussed by 
M. Znojil jl9| I2IH l2T| and A. Mostafazadeh [22j. From a phenomenological point of 
view, I. Rotter discussed double poles in the scattering matrix and level repulsion of 
unbound states in the framework of an effective many body theory fl^ 1^ |23 12E] • 

The characterization of the singularities of the eigenenergy surfaces at a degeneracy 
of unbound states in parameter space arises naturally in connection with the Berry phase 
of unbound states that was predicted by Hernandez et al |2Z1 EHl I2nj and independently 
by W.D. Heiss jSO]; see also the interesting recent work by A. A. Mailybaev et al [HI]. The 
Berry phase of two unbound states was measured by the Darmstadt group jHSESI- The 
unfolding of energy eigenvalue surfaces at a degeneracy of unbound states of a Hermitian 
Hamiltonian was discussed by E. Hernandez et al [33] and the unfolding of eigenvalue 
surfaces of non-Hermitian Hamiltonian matrices with applications in modern problems 
of quantum physics, cristal optics, physical chemistry, acoustics, mechanics and circuit 
theory has been the subject of many recent investigacions (35 ] 136 ] l37 | l38 | l39 t lln | HT ] Il2] . 

In this paper, we will be concerned with some physical and mathematical aspects 
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Figure 1. The double barrier potential is such that it supports two unbound states 
with equal energies and half- lives. The control parameters of the system are d and V3. 



of the mixing and degeneracy of two unbound energy eigenstates in an isolated doublet 
of resonances of a quantum system depending on two control parameters. The plan of 
the paper is as follows: In section 2, we give the results of a numerical computation of 
the surfaces that represent the resonance energy eigenvalues as functions of the control 
parameters in the scattering of a beam of particles by a double barrier potential. The 
analytical structure of the singularity of the energy surface at the crossing point is 
characterized in section 3, where we also introduce a contact equivalent approximant 
to the energy surface at the degeneracy point. Section 4 is devoted to a discussion of 
crossings and anticrossings of energies and widths, as well as, the changes of identity of 
the poles of the jS"— matrix, in term of sections of the energy surfaces. We end our paper 
with a short summary and some conclusions. 

2. Resonances in a double barrier potential 

Doublets of resonances and accidental degeneracies of unbound states may occur in 
the scattering of a beam of particles by a potential with two regions of trapping. A 
simple example is provided by a spherically symmetric potential V{r) such that the two 
regions of trapping are two potential wells defined by two concentric potential barriers 
located between the origin of coordinates and the outer region where the potential V{r) 
vanishes. In order to make the analysis as simple and explicit as possible, we take the 
wells and barriers to be square as shown in figure 1. In this section, we will consider the 
conditions for the occurrence of a degeneracy of unbound states in this simple system 
and we will describe the numerically exact computation of the surfaces that represent 
the complex energy eigenvalues as functions of the control parameters of the system in 
the neighbourhood of and at a degeneracy of unbound states. 
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2.1. The Jost regular solution 

The s-wave radial Schrodinger equation is 



d^(j){k, r) 
dr"^ 



+ - U{r))(P{k,r) = 0, 



(1) 



the potential U{r) = 2mV{r)/fi^ is a double barrier such that between the origin of 
coordinates and the outer region, where the particles propagate freely, there are two 
square potential wells separated by two square potential barriers, as shown in figure 
1. The system has seven parameters, the positions rj (i = 1,2,3,4) and heights Vi (i 
= 2,3,4) of the four discontinuities of the potential. We will keep the five parameters 
(V2, V4,ri,r3 — r2,r4 — r^) fixed and will vary the depth V3 of the outer well and the 
thickness of the inner barrier d = r2 — ri. In the following, we will refer to the pair of 
parameters (rf, V3) as the control parameters of the system. 

The radial Schrodinger equation ((T)) is solved exactly. The Jost regular solution of 
normalized to unit slope at the origen, (j){k,r), is as follows: 

In the wells, 

01 (/c, r) = — sin /cr, < r < ri, (2) 



k 



and 



4>3{k,r) = 4>2{k,r2) cos(^K3{k,V3){r - r2 
+ a2{k,d) sm(K3{k,V3){r - r2) 



r2 < r < r^. 



(3) 



In the barriers. 



(t)i{k,r) = 0i_i(fc,ri_i) cos\i(^Ki{k){r - ri^i)^ 



+ at 



.i{k,d) smh( Ki{k){r - 



n-i, 



Ti-i <r <ri 



2,4, 



and, in the outer region, 

05(fc, r) = 04(/c, r^) cos k{r — r^) + a^^ik, d) sin k{r — r^) 

In these expressions k is the wave number of the free waves and 

1/2 



K,[k) = ({-i)\Ui-e)) 



2,3,4 



(4) 

r4 < r < 00. 

(5) 

(6) 



is the wave number in the barriers and the outer well. 

The logaritmic derivatives ai(/c,(i, V3) of (f>{k,r) at the consecutive discontinuities 
ri and rj+i are related by the matching conditions at rj+i, 

k 



ai{k) 



K2{k) a^{k)+timh{K2{k)d) 
■cotKri, a2{k, d) = ttt , n y'l 



a3{k;d,V3 



K2{k) ' ^ ' ' K3{k)l + ai{k)ia.Yi\i{K2{k)d) 

Ksik, Vs) a2{k, d) - tan{K,{k){rs - rs)) 



Ki{k) 1 + a2{k, d) tan{K3{k){r3 - rs)) ^ 



(8) 
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and 



a4(/c; d, 1/3) 



K^{k) a^{k- d, V3) + tanh(i^4(A;)(r4 - r^)) 



(9) 



k 1 + a3(fc; d, V3) tanh(_ft'4(A;)(r4 — ra)) 

Since the first logaritmic derivative is known, succesive substitution of ai{k,d,Vs) 
in ai+i{k, d, V3) gives an explicit solution for a^^k; d, V3) and an explicit expression for 
the regular solution is obtained from equations (j21ISI)- 

2.2. The Jost function 

When the regular solution in the outer region, is written as a linear combination of 
an outgoing wave exp{ikr) and an incoming wave exp {—ikr) 

05(fc,r) = 04(/c,r4)^ ^1 - iaiik; V3,d)^ expik{r - r^) 
+ ( 1 + iai{k] V3, d) ) exp —ik{r — 



T'i < r < 00, 



(10) 



the coefficient of the incoming wave is the Jost function. Making use of equations (0}, 
Q and (fTUj) we find an exact expression for the Jost function f{—k), 

f{—k]d,Vz) = sinA;ri cos\iK2{k)d + ai{k) smh.K2{k)d 



X 



—I 



cos Ksi^k, V3){r3 - ra) + a2{k, d, V3) sin K3{k, V3){r3 - ra) 
K,{k) 



k 



sinh _ft'4(fc)(r4 — r^) + 0:3 (fc; d, V3) cosh_ft'4(A;)(r4 — r^) 



coshi^r4(A;)(r4 — r^) + a3(/c; d, V3) sinhK4(/c)(r4 — r3)| | expiA;r4 (11) 



2.3. The physical solutions 

The scattering wave function, ip{k,r), and the regular solution (f){k,r) are related by 



ipi^k, r) 



-2ik 



(f){k,r) 



(12) 



(13) 



fi-k) 

and the scattering matrix is given by 

where the Jost function f{—k) is given by (fTT|) . 

The zeros of the Jost function give resonance poles in the scattering wave function 
ip{k, r), and in the S{k) matrix, and from (Q) and (j21EI); we may verify that all unbound 
energy eigenfunctions of the radial Schrodinger equation are associated with roots (zeros) 
of the Jost function. 

Unbound state eigenfunctions also called resonant-state or Gamow eigenfunctions 
are the solutions of equation that vanish at the origin, and at infinity satisfy the 
outgoing wave boundary condition. When the Jost function has 3j zero dbt , th-G 
coefScient of the incoming wave in (ITUj) vanishes, and (p{kn,r) is proportional to the 
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outgoing wave solution of equation (Q) for r larger than the range of the potential. 
Hence, the unbound state eigenf unctions are related to the regular solution by 

Un{kn,r)=N-^(i){kn,r), (14) 

where is a normahzation constant and /c„ is a zero of the Jost function 

/(-A;„;d,\/3) =0 (15) 

2.4- Degeneracy of unbound energy eigenvalues 

A degeneracy of unbound states results from the exact coincidence of two simple zeros 
of the Jost function, which merge into one double zero lying in the fourth quadrant of 
the complex fc— plane. Hence, the condition for a degeneracy of two unbound energy 
eigenstates at some k = kd is that both, the Jost function and its first derivative vanish 

at kj,, 

where f{—k; d, V3) is given in (fTTj). 

Therefore, to locate a degeneracy of unbound states, the coupled equations (fTT)|) 
were solved numerically. The zeroes of the Jost function are found by an algebraic 
computer package that searches for the minima of \ f{—k)\ in the complex k-plane. In 
the computation the five parameters V2, V4,ri,r3 — r2 and r4 — r^ were kept fixed at 
the values V2 = V4 = 2, ri = 1, — r2 = l,r4 — r^ = 0.304892 and only the control 
parameters d and V3 were allowed to vary. Starting with the values d = 2 and V3 = 

1.04, we find the first isolated doublet of resonances at ki = 2.2101546 - i 0.1366887 and 
k2 = 2.2321776 - i 0.0017984. 

Then, we adjusted the control parameters d and V3 until ki and k2 became equal 
to some common value k^- We also computed numerically \df{—k)/dk\ aX k = k^ 
to verify that the second equation was also satisfied to some previously prescribed 
accuracy. In this way, we found that by fine tuning the control parameters to the 
values d* = 1.1314661145 and = 1.038235081, the first doublet of resonances 
becomes degenerate, with a precision better than one part in 10^, at k^ = 2.22697606 - 
i 0.07220139. 

In the following we will refer to the double zero of f{—k) at k^ as the degeneracy 
point or crossing point of the doublet of resonant states in the complex k-plane and to 
the point {d*, V^) as the exceptional point in parameter space. 

2.5. Energy surfaces 

The energy eigenvalues £n{d, V3) = {h'^ /2m)k'^{d, V3) of the physical system are obtained 
from the zeroes of the Jost function, given in equation (fTH|) . That condition defines, 
implicitly, the inverse functions 

kn{d, V3) = f-\0; d, V3), n = 1, 2, .... (17) 
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Figure 2. The two - sheeted surface Sr that represents the real part of the eigenwave 
numbers ki and ^2 as functions of the control parameters (d, V3) in the vecinity of a 
degeneracy of unbound states. 

as branches of a multivalued function ^3] which will be called the wave number pole 
position function. Each branch kn{d, V3) is a continuous, single valued function of the 
control parameters. When the physical system has an isolated doublet of unbound 
states which become degenerate for some exceptional values of the control parameters 
(d*, VJj*), the corresponding two branches of the energy pole position function, £i{d, V3) 
and £2{,d, V3), are equal (cross or coincide) at that point. With the purpose of exploring 
the geometrical and topological properties of the surfaces representing the pole position 
function fci,2(c?, V3) of the isolated doublet of resonances in parameter space, we solved 
numerically the implicit equation p5|) for ki{d,V^) and /c2((i, V3) in the neighbourhood 
of and at a degeneracy of unbound states. The results of the numerical computation 
are represented as surfaces in a Euclidean space with coordinates (Re A;, Im /c, d, V3). 

In figure 2, the real function Re ki^2{d, V3) is shown as a surface Sr in the three- 
dimensional subspace with cartesian coordinates (Re k, d, V3). Similarly, in figure 3, the 
real function Im ki 2{d, V3) is shown as a surface Sj in the three-dimensional subspace 
with cartesian coordinates (Im k, d, V3). 

We see that close to the degeneracy of unbound states, the function Re /ci,2('^, V3) 
has two branches and is represented by a two sheeted surface Sr. The two sheets of Sr 
are two copies of the plane {d, V3) which are cut and joined smoothly along a line Cr 
starting at the exceptional point and extending to values d > d* and V3 > I3*. 

The function Im ki^2{d, V3) has also two branches and is represented by a two sheeted 
surface Sj. The two sheets of 5"/ are two copies of the plane {d, V3) which are cut and 
pasted smoothly along a line Cj extending from the exceptional point to values d < d* 
and V3 < V^. 

The projection of the lines Cr and Cj on the plane [d, V3) are the two halves of a 



Non-Hermitian degeneracy of two unbound states 



8 




Figure 3. The two-sheeted surface Si that represents the imaginary part of eigenwave 
numbers fci and fc2 as functions of the control parameters (d, V3) in the neighbourhood 
of a degeneracy of unbound states with complex resonance energies Ei = h?kf/2m, 
i = 1,2. 

line C that goes through the exceptional point {d*, V^). 

At the exceptional point, and only at that point, both the real and imaginary parts 
of ki and k2 are equal. Therefore, in the complex k-plane, at the crossing point, the two 
simple zeroes of the Jost function merge into one double zero which is an isolated point 
in the complex k-plane. 

3. The analytical behaviour of the pole-position function close to the 
exceptional point 

In this section, it will be shown that, the singularity of the energy surface at the 
degeneracy or crossing point is an algebraic branch point of square root type. The 
analytical structure of the singularity of the energy surface will be worked out and 
discussed in detail. The results in this and the following section are not restricted to the 
double barrier potential, to stress this generality, the control parameters will be called 
Xi and X2. 

Let us start by recalling that, when the first and second absolute moments of 
the potential exist, and the potential decreases at infinity faster than any exponential 
[e.g., if V{r) has a Gaussian tail or if it vanishes identically beyond a finite radius], 
the Jost function f{—k;xi,X2), is an entire function of k |43j. The entire function 
of k, /(— /c; Xi, X2), may be written in the form of an infinite product, according to 
Hadamard's form of the Weierstrass factorization theorem ^1], and by using a theorem 
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of Pfluger ; we may write 

k 



f{-k; Xi,X2) = /(0)exp(i/ci?) TT( 1 - — 



R is the range of the potential, /(O) 7^ and are the zeroes of f{—k] xi,X2), see 
also R.G. Newton gg. 

When the Jost function has a set of isolated zeroes, the implicit equation for the 
pole position function, equation (|TKjl . may, in principle, be solved, at least numerically, 
for each branch kn{xi,X2) without ambiguity. When the system has an isolated doublet 
of resonances which may become degenerate, the corresponding branches of the pole 
position function, say ki{xi,X2) and ^2(3^1, 3^2), may be equal (cross or coincide) at an 
exceptional point. In this case, it is not possible to solve equation (|THjl for each individual 
branch without ambiguity and one should proceed to solve that equation for the pole 
position function of the two members of the isolated doublet of resonances. 

To be precise, we will say that the system has an isolated doublet of unbound 
states if there is a finite bounded and connected region Ai in parameter space and 
a finite domain V in the fourth quadrant of the complex A;— plane, such that when 
(xi,X2) e Ai, the Jost function has two and only two zeroes, ki and k2, in the finite 
domain V C C, all other zeroes of f{—k] xi,X2) lying outside V. 

The pole position function ki^2{xi,X2) of the isolated doublet of resonances, 

^1,2(3:1, X2) = /~^(0; xi, X2), (19) 
is implicitly defined by the equation f{—ki^2',xi,X2) = 0, and the conditions 

df{-k;x„X2U ^d^f{-k,x„X2h (20) 



dk J ka ^ V dk"^ J ka 

for {xi,X2) in a neighbourhood of the exceptional point (x^,X2). 

Then, to find an expression for the pole position function of the isolated doublet of 
unbound states, the two zeroes of /(— fc; xi, X2), corresponding to the isolated doublet 
of unbound states are explicitly factorized in (fTH|) as 

f{-k;xi,X2) = k~ki{xi,X2) k - k2{xi, X2) gi,2ik; Xi, X2), (21) 

which may be conveniently rearranged as 

f{-k; Xi,X2) = ik-^{ki + k2) f -liki-k2f gi,2ik,Xi,X2). (22) 

with 

gi,2ik; Xi, X2) = /(O) exp{ikR) J] " j, (r r ) ) ' 

Ki[Xi, X2)K2[Xi, X2) g V Kn[Xi,X2)^ 

(23) 

the expression in square brackets in the right hand side of equation ()22|) is the Weierstrass 
polynomial of the isolated doublet of unbound states [iH]. 

Solving equation for fci,2(a;i, ^2) when f{—k;xi,X2) vanishes, we get 
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kl,2{Xl, X2) = ^ (^h{Xl,X2) + k2{Xi,X2)^ + )J^(^ki{xi,X2) - k2{Xi,X2)^ 

(24) 

with {xi,X2) in a neighbourhood of the exceptional point. This equation relates the 
wave number pole position function of the doublet to the pole position function of the 
individual unbound (resonance) states. Since the argument of the square-root function 
is complex, it is necessary to specify the branch. Here and thereafter, the square root 
of any complex quantity F will be defined by 



F\ ex.p{i-argF), 



< arg F < 27r 



(25) 



so that IV F\ = ^^/\F\ and the F— plane is cut along the real axis. 

Now, we will proceed to the derivation of a contact equivalent approximant to the 
pole position function of the doublet at the crossing point. 

According to the preparation theorem of Weierstrass jlH], the functions 
1/2(^1 (xi, 0:2) + k2{xi, X2)) and l/4:{ki{xi,X2) — k2{xi, ^2))^, appearing in the right hand 
side of equation ()24|) . are regular at the exceptional point and admit a Taylor series 
expansion about that point. 



^(^ki{xi,X2) + k 



■2[Xi,X2) ) = kd{xl,x*2) 
2 



i=l 



and 



{ki{xi,X2) - k2{xi,X2)^ = ELi4^^(a^i - a;*) + 0((Xi -x*)2). 



(26) 



(27) 



The complex coefficients cf^ and df\ appearing in these equations, may readily be 
computed from the Jost function with the help of the implicit function theorem 



.(1) 



df(—k;x\,X2) 
dxi 



X2-i k=kH. 



(28) 



and 



d] 



(1) 



f{-k:xi,X2) 



f d^f{-k;xi,X2) \ 
\ dxidk J 



X2' 



k=kj 



S'^f(-k;xi,x2) 



X 



2 

df{-k;xi,X2) 
dxi 



f{—k\xi,X2) 



k=kj 



(29) 



X2' 



A contact equivalent approximant, fci, 2(2^1, ^2), to the doublet's pole position function 
is obtained when the Taylor series expansions and (|77|) are substituted for the 
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functions l/2(A;i(xi, X2) + k2{xi,X2)) and l/4:{ki{xi,X2) — ^2(2^1) 3^2))^- Keeping only the 
first order terms, we obtain 

2 

i=l 

+ ^J\[c?ixi-xl) + c^^\x2-x*2)]. (30) 

Tfien, from (j221) and (jSDI), close to the exceptional point, the Jost function may be 
approximated as 

f{-k; xi, X2) ^ 77-—* — T^fdoubi-k] xi, X2), (31) 

gi,2[kd]x\,xl) 

where 

1 2 



fdoub{-k; xi,X2) = k - (kd + y^^d^l\xi,x*)) 

i=l 



-l{c?ix,-xl) + c^\x2-x;)). (32) 

The coefficient [(71,2(^^5 a^i, 2^2)]"^ "^^Y be understood as a scaling factor. 

Hence, the two parameter family of functions fi,2{~k',Xi,X2) is contact equivalent to 
the Jost function at the exceptional point and is also a universal unfolding [47j of 
/(— fc; Xi,X2) at the exceptional point where the degeneracy of unbound states occurs. 

3.1. Energy-pole position function 

A contact equivalent approximant to the energy pole position function <£^i, 2(2^1, 2^2) at 
the crossing point of the doublet of unbound states is readily obtained from ()24p2(jl - 
fI7\i . Taking the square in both sides of (jHUj) . multiplying them by /2m and recalling 
£i = {h?/2m)kf, in the approximation of ()2(ill3()j) . we get 

^^1,2(3:1, X2) ^ £d{xl, xl) + A£d{xi, X2) + 61,2(3^1, X2), (33) 

where 



ei,2(xi,X2) = )J\[Ci'\x, - xt) + Ci'\x2 - x*2)], (34) 
and Cl^'^ = (h'^kd{xl,X2)/n2) cf'K It will be convenient to change slightly the notation 



The components of the real fixed vectors R and / are the real and imaginary 
parts of the coefficients C^^^ of (x, — x*) in the Taylor expansion of the function 
l/A{Si{xi,X2) — £^2(3^1,3^2))^ and the real vector ^ is the position vector of the point 
(xi,X2) relative to the exceptional point (x*,X2) in parameter space. 
In the notation defined in equations (jHSj) . 

elJx,,X2) = l(iR-0+^iI-0) (36) 
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and 



Solving for the real and imaginary parts of the function ei,2(a;i, X2), we obtain 

1/2 



Re £1,2(2:1, X2) 
Im eio{xi,X2) 



"2V2 
1 

"271 



(i?-e? + (/-o'-i?-e 



.11/2 



and 



sgn{^Reei^2^sgn{^mti^2^ = sgn(^I-^. 



(37) 

(38) 
(39) 

(40) 



It follows from (j38jl . that Re ei,2(a;i, 0:2) is a two branched function of (G; ^2) which may 
be represented as a two-sheeted surface Sr in a three dimensional Euclidean space with 
cartesian coordinates {Re£,^i,^2)- The two branches of Re£{C,i,^2) are represented by 
two sheets which are copies of the plane (^1, ^2) cut along a line where the two branches 
of the function are joined smoothly. Since a negative and a positive numbers are equal 
only when both vanish, the cut is defined as the locus of the points where the argument 
of the square- root function in the right hand side of (jHHj) vanishes. Close to the origen 
of coordinates (the exceptional point), this locus is defined by a unit vector in the 

— * — > 

(^15^2), plane such that 

/ ■ io = and R-io = -\R- L\ (41) 

Therefore, the real part of the energy-pole position function, Si,2{xi,X2), as a 
function of the real parameters (a;i,a;2), has an algebraic branch point of square root 
type (rank one) at the exceptional point with coordinates (x^,X2) in parameter space, 
and a branch cut along a line, Cn, that starts at the exceptional point and extends in 
the positive direction defined by the unit vector satisfying equations (jH]). 

A similar analysis shows that, the imaginary part of the energy-pole position 
function, Im i^i,2(a;i, ^2), as a function of the real parameters (a;i,a;2), also has an 
algebraic branch point of square root type (rank one) at the exceptional point with 
coordinates (x^, Xg) in parameter space, and also has a branch cut along a line, £/, that 
starts at the exceptional point and extends in the negative direction defined by the unit 
vector satisfying equations (jH}. 

The branch cut lines, Cr and Cj, are in orthogonal subspaces of a four dimensional 
Euclidean space with coordinates {ReS, ImS, ^i, ^2) - but have one point in common, 
the exceptional point with coordinates ( 

Along the line £r, excluding the exceptional point (x*,X2), Re Si{xi,X2) = Re 
£^2(3^1, X2) but Im £^1(0:1,0:2) 7^ Im £2{xi,X2). 

Similarly, along the line Cj, excluding the exceptional point, Im £^i(xi,X2) = Im 
£2{xi,X2), but Re £^i(a;i,a;2) 7^ Re £2{xi,X2). 
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Figure 4. The curves C'„(7ri) and C'„+i(7ri) are the intersection of the hyperplane 
TTi : ^2 = £,2^ and the two-sheeted surface tn,n+i- The projections of Cn{T^i) 
and C'n+i(7ri) on the planes {Im£,£i) and (Re5,^i) show a crossing of widths and 
anticrossing of energies, respectively. The projections of C„(7ri) and C„+i(7ri) on the 
plane {Re£, Im£) are the trajectories of the 5— matrix poles in the complex energy 
plane. In the figure, = d — d* . 

Equality of the complex resonance energy eigenvalues (degeneracy of resonances), 
£i{xl,X2) = £^2(3^1,3^2) = ^d{xl,X2), occurs only at the exceptional point with 
coordinates (x^,X2) in parameter space and only at that point. 

In consequence, in the complex energy plane, the crossing point of two simple 
resonance poles of the scattering matrix is an isolated point where the scattering matrix 
has one double resonance pole. 

4. Phenomenology of the exceptional point 

There is a variety of phenomenological manifestations of the topological and geometrical 
properties of the singularity of the energy surfaces at an exceptional point, which fall 
roughly in two categories, according to the experimental set up. First, when one control 
parameter is slowly varied while keeping the other constant, cossings and anticrossings 
of energies and widths are experimentally observed jlHl IIHI EHI , as well as, the so called, 
changes of identity of the poles of the scattering matrix [HIEI]. Second, when the system 
is slowly transported in a double circuit around the exceptional point in parameter space, 
it is observed that the wave function of the system acquires a geometrical or Berry phase 

1211 EHl EHl EDI EB Ea ES! . 
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Figure 5. The curves C'„(7r2) and Cn+i{T^2) are the intersections of the hyperplane 
1^2 that goes through the exceptional point (^i,^!) parameter space and the two- 
sheeted surface ei,2('^i, ^2)- The projections of Cn{'^2) and Cn+i{i^2) on the planes 
{Re£,^i) and {Im£,^i) show a joint crossing of energies and widths. The projections 
of C'„(7r2) and C'„+i(7r2) on the plane {Re£,Im£) are the critical trajectories of the 
5— matrix poles in the complex energy plane. At the crossing point, the two simple 
poles coalesce into one double pole of S{E). 



4-1- Sections of the energy surfaces 

The experimentally determined dependence of the difference of complex resonance 
energy eigenvalues on one control parameter, say ^1, while the other is kept constant, 
5 = #,i = 1,2,3, 

^1(6, eT) - 4(6, = 61,2(6, g'^) (42) 
has a simple and straightforward geometrical interpretation, it is a direct measurement of 
the intersection of the eigenenergy surface of the doublet £1,2(6, 6) with the hyperplane 
defined by the condition ^2 = g"*, i=l,2,3. 

The intersection of the eigenenergy surface €1,2(6,6) ^^'^ each one of the 
hyperplanes TCi ■ ^2 = g^ , defines two three-dimensional curves for each value of g^ 

f C'i(7ri) 

M6,6)n7r,= (43) 
(^2(,7rjj 

The sections Ci^Hi) and 6*2 (vTj) are the three-dimensional curves traced by the points 
£1(6, g'*) and 4(6, g"*) on the surface 61,2(6,6) when the point with coordinates 
(6,^"*) moves along a straigth line path parallel to the O^i axis, and ^i,j < 6 ^ 6,/, 
in parameter space. 
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Figure 6. The curves Cni'^'i) and C'„+i(7r3) are the mtersection of the hiperplane 
TTa := ^2 and the two sheeted surface e„,„+i. The projections of C„(7r3) and C„+i(7r3) 
on the planes (Re ^i) and (Im £, ^i) show a crossing of energies and an anticrossing 
of widths. The projections of C„(7r3) and C„+i(7r3) on the plane {Re£,Im£) are the 
trajectories of the 5— matrix poles in the complex energy plane and they do not cross. 
In the figure, = d — d* . 



The projections of the sections Ci{T^i) and C2{'^i) on the planes (Re and 
(Im^, ^i) are 

RelCmi-n,)] = ReSmi^,, g*)) m = 1, 2 (44) 

and 

ImlC^n,)] = ImEmiii, m = 1, 2 (45) 

respectively, see figures 4-6. A comparison of the representations of the eigenenergy 
surfaces provided by the numerically exact computation and the contact equivalent 
approximant is shown in figure 7. 

The projections of the sections CiijTi) and C2{T^i) on the plane {Re£, Im £) are the 
trajectories of the 5— matrix poles in the complex energy plane. An equation for these 
trajectories is obtained by eliminating between Re^m(^i, ^2*'') Im^m(^i, .^2*'')) ^ = 
1,2, eauations dHHm . 

(Re£m) ' - 2 cot 01 (Re8^) [imSj^ - [imSj^ ^ -\[R- Ic'^) = (46) 



where 



cot0i = ^ (47) 
-'1 



and the constant vector is such that 
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Figure 7. Projections of the sections Cn{T^i) and C'„+i(7ri) on the planes (Im 
and (Re £,S,i) are shown in the left and right columns respectively for a) V3 — 1.0381, 
b) V3 = V3* and c) V3 =1.0834. The ftill line is the numerically exact calculation, the 
dotted line is the contact equivalent approximant. 



The discriminant of P6|l . 4 cot^ 0i + 1, is positive. Therefore, close to the crossing point, 
the trajectories of the S*— matrix poles are the branches of a hyperbola defined by 

4-2. Crossings and anticrossings of energies and widths 

Crossings or anticrossings of energies and widths are experimentally observed when the 
difference of complex resonance energy eigenvalues, 'C2)~^2('Ci) ^2) = A^E — zl/2Ar 

is measured as a function of the slowly varying parameter ,^1, keeping the other constant. 
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f(0 

S2 ■ 



From equations 

AE = Er. - 



, and keeping ^2 = ^2*^ obtain 



^n+l — — 



and 



AF 



F„ -F 



n+l 



-V2 



6= 

1/2 



(49) 



■(») 



(50) 



These expressions allow us to relate the terms (J-^) directly with observables of 

the isolated doublet of resonances. Taking the product of AEAT, and recalling equation 

we get 



AEAT 



■(i) 



(51) 



and taking the differences of the squares of the left hand sides of and we get 



(52) 



At a crossing of energies AE vanishes, and at a crossing of widths AF vanishes. 
Hence, the relation found in eq. (j51|) means that a crossing of energies or widths can 
occur if and only if (/ ■ Oc(') vanishes 

For a vanishing (/•^c)f(o = = Ai?AF, we find three cases, which are distinguished 

by the sign of (R ■ ^c)c(»)- From eqs. (11^ and (1301) . 

(i) (R ■ ^c)t(i) > implies AE ^ and AF = 0, that is, energy anticrossing and width 

S2 

crossing. 

(ii) (i?-^c)c(0 = implies AE = ano? AF = 0, that is, joint energy and width crossings, 

?2 

which is also degeneracy of the two complex resonance energy eigenvalues. 

(iii) (R ■ ^c)c(i) < implies AE = and AT 7^ 0, i.e. energy crossing and width 

S2 

anticrossing. 

This rich physical scenario of crossings and anticrossings for the energies and widths 
of the complex resonance energy eigenvalues extends a theorem of von Neumann and 
Wigner [5] for bound states to the case of unbound states. 

The general character of the crossing-anticrossing relations of the energies and 
widths of a mixing isolated doublet of resonances, discussed above, has been 
experimentally established by P. von Brentano and his collaborators in a series of 
beautiful experiments |IH t IHH]. 



4-3. Trajectories of the S— matrix poles and changes of identity 

In subsection 4.1, we found that, when one control parameter, say ^1, is varied and 
the other control parameter is kept constant ^2*'' close to the exceptional value, the 
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trajectories of the S"— matrix poles are the branches of a hyperbola defined by 
The asymptotes of this hyperbola are the two straight lines defined by 



= tan ^Re£^^^ (53) 



and 



= _ cot |i?e^(^^) (54) 

The two branches of the hyperbola are in opposite quadrants of the complex energy 
plane divided by the asymptotes, see figure 8. 

We find three types of trajectories, which are distinguished by the sign of (i?-^) |^2=6 • 

(i) When {R ■ Ol6=6 ^ 0' branch of the hyperbola lies to the left and the other 
lies to the right of a vertical line that goes through the crossing point. 

(ii) Critical trajectories, when {R ■ Ol6=6 ~ 0' trajectories are the asymptotes of 
the hyperbola. The two simple poles start from opposite ends of the same straight 
line and move towards each other until they meet at the crossing point where they 
coalesce to form a double pole of the S"— matrix. From there, they separate moving 
away from each other on a straight line at 90° with respect to the first asymptote. 

(iii) When {R ■ Ol6=6 ^ '^^^ P°^^ moves on one branch of the hyperbola that hes 
above and the other pole moves on the other branch that lies below a horizontal 
straight line that goes through the crossing point. 

When a small change in the control parameter ^2*'' changes the sign of {R ■ 016=6' 
it produces a small change in the initial position of the poles, but the trajectories 
change suddenly from type (i) to type (iii), this very large and sudden change of the 
trajectories exchanges almost exactly the final position of the poles as can be appreciated 
from figure 8. This dramatic change was called a "change of identity" by Vanroose, van 
Leuven, Arickx and Broeckhove |H| who discussed an example of this phenomenon in the 
5— matrix poles in a two-channel model, Vanroose has also discussed these properties 
in the case of the scattering of a beam of particles by a double barrier potential with 
two regions of trapping. 



4.4- Changes of identity when going around the exceptional point 

In the previous discussion of the trajectories of the S"— matrix poles and their changes 
of identity, the two straight lines, parallel to the O^i axis, defined by the conditions 
< ^1 < ^ij and ^2 = ^2*^ i = 1;3, are the two long sides of a very long and narrow 
rectangle in parameter space that surrounds the exceptional point. The two short sides 
of this rectangle are defined by the small change in ,^2 when going from ^2^^ to 
keeping or ^ij constant, to which we refered above, when explaining the changes of 
identity of the poles of the doublet of unbound states. In order to understand better 
the meaning of the changes of identity of those poles, it will be convenient to deform 
continuously the long and narrow rectangle into a circle and consider the motion of the 
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Figure 8. Trajectories of the poles of tlie scattering matrix, S{£), of an isolated 
doublet of resonances in a double barrier potential close to a degeneracy of unbound 

states. The control parameters are the width, rf, of the inner barrier and the depth, 
V3, of the outer well. The trajectories are traced by the poles f 1 (d, VJj*'*' ) and 
£2{d,V^^^) on the complex f— plane when the point {d,V^^^) moves on the straight 
line path tt^; V3 = V^^\ The top, middle, and bottom figures show the trajectories 
corresponding to [R ■ Cc),(i) > 0, {R ■ Cc)t(i) = 0, and [R ■ ^c)f{i) < 0, respectively, 
with (^1, ^2) = {d^ d* ,V^~V^). The full line is the numerically exact calculation, the 
dotted line is the contact equivalent approximant. 
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Figure 9. The two poles of the isolated doublet of unbound states trace a starlike 
trajectory in the complex fc— plane, shown on the left hand side of the figure, when the 
physical system is transported in parameter space twice around the exceptional point 
in the circular path shown at the right hand side of the figure. For an explanation, see 
the text. 

zeroes of the Jost function A;i(xi, ^2) and k'2,{xi,X2) in the complex /c— plane, when the 
system is transported in parameter space around the exceptional point in the closed 
circular path equivalent by continuous deformation to the long and narrow rectangle of 
the previous discussion. 

On the left hand side of figure 9, we show the trajectory traced in the complex 
A;— plane by the two zeroes of the Jost function, ki{xi^X2) and k2{xi,X2)^ when the 
physical system is transported in parameter space twice around the exceptional point 
on the double circular circuit shown on the right hand side of the same figure. When 
the system is at the crossing point of the two circular paths in parameter space, one 
zero, say ki{xi,X2)i is at the uppermost point and the other, say A;2(a^i, 2^2), is at the 
lowermost point of the four points star in the complex /c— plane. 

(i) As the point representing the system in parameter space moves on the inner circle 
in counterclockwise direction, from its initial position until it makes one complete 
round about the exceptional point and is back at its initial position, the ki{x\, X2) 
zero moves on the star in the complex A;— plane in clockwise direction, from its 
initial position at the topmost point of the star, goes through the point at the 
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extreme right hand side of the star and ends at the lowest point in the star, while 
the ^2(3^1, X2) zero moves also in clockwise direction, from its initial position at the 
lowest point in the star goes trough the point in the extreme left hand side of the 
star and ends at topmost point on the star. It follows that, when the system goes 
around the exceptional point once in parameter space, the poles of the scattering 
matrix are exactly exchanged. 

(ii) As the point representing the system in parameter space continues its 
conterclockwise motion, now on the outercircle, until it completes a second round 
about the exceptional point and is back at its initial position, the A;i(a;i,a;2) zero 
moves on the star in the complex /c— plane in clockwise direction from the lowest 
point on the star, goes through the extreme left hand side point on the star and 
ends at the topmost point on the star, while the k2{xi,X2) zero moves on the star 
in the complex fc— plane also in clockwise direction, from the topmost point on 
the star, goes through the extreme right hand side point on the star and ends at 
lowest point of the four points star. Therefore, when the system goes around the 
exceptional point twice in parameter space the poles of the S{k) matrix and the 
complex energy eigenvalues return to their initial values in the complex /c— plane. 
The eigenfunctions also return to their initial values but they acquire a geometric 

phase isziEHiEniisniisiiES. 

5. Summary and conclusion 

In this paper, we discussed some mathematical and physical aspects of the non- 
Hermitian degeneracy of two unbound energy eigenstates of a Hamiltonian dependending 
on two control parameters. We solved numerically the implicit transcendental equation 
that defines the eigenenergy surface of a degenerating isolated doublet of unbound states 
in the simple but illustrative case of the scattering of a beam of particles by a double 
square barrier potential. The analytical characterization of the singularity of the energy 
surface was made in the more general case of a short ranged potential with two regions of 
trapping. We showed that, from the explicit knowledge of the Jost function as a function 
of the control parameters of the system, it is possible to derive a two parameter family 
of functions which is contact equivalent to the exact energy-pole position function at 
the degeneracy point and includes all small perturbations of the degeneracy conditions. 
This unfolding of the degeneracy point gives a simple and explicit, but very accurate, 
representation of the eigenenergy surface close to the exceptional point, see figure 7. 
In parameter space, the surface that represents the complex energy eigenvalues has a 
branch point of square root type at the crossing point, and branch cuts in its real and 
imaginary parts that start at the exceptional point but extend in opposite directions 
in parameter space. In the complex energy plane, the crossing point of two simple 
resonance poles of the scattering matrix is an isolated point where the scattering matrix 
has one double resonance pole. Crossings and anticrossings of the energies and widths 
of the resonances in an isolated doublet of unbound states of a quantum system, as well 
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as the sudden change in the shape of the 5"— matrix pole trajectories, observed when 
one control parameter is varied while the other is kept constant at a value close to the 
exceptional value, are fully explained in terms of sections of the energy surfaces. 
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